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Abstract 

The superposition of two independent point processes can be described by 
multiplication of their probability generating functionals (p.g.fi.s). The inverse 
operation, which can be viewed as a deconvolution, is defined by dividing the 
superposed process by one of its constituent p.g.fi.s. The deconvolved process 
is computed using the higher-order chain rule for Gateaux differentials. The 
higher-order quotient rule for Gateaux differentials is first established and then 
applied to point processes. 
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The superposition of two point processes can be formed by multiplication of their 
probabihty generating functionals [5] • While it is clear that the inverse operation would 
involve division of the probability generating functionals, there are no results regarding 
the Janossy densities after this operation. This is probably due to the fact that the 
denominator involves an infinite sum and there was no easy way of computing the 
point process. In this article, I derive the deconvolved point process using the recently 
derived higher-order chain rule for Gateaux differentials jT|. 

The paper is structured as follows. Section 1 describes the necessary background on 
functionals to derive the result, including the probability generating functional (p.g.fl.), 
the Gateaux differential, Janossy densities and the higher-order chain rule. Section 2 
presents the main result by first establishing the n*''-order Gateaux differential of the 
reciprocal of a functional in Lemma 4. Theorem 1 establishes the higher-order quotient 
rule for Gateaux differentials using Leibniz' rule and the recently derived Faa di Bruno's 
formula for Gateaux differentials [T]. This is used to determine the expression for the 
deconvolved point process via the probability generating functional [6 . 

1. Probability generating functionals 

This section discusses generating functionals as a means of describing systems with 
a variable number of objects. Gateaux differentials are discussed as a means of 
detemining their constituent functions. The higher-order chain rule is introduced to 
determine the differentials of composite functionals [T] . 
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Definition: probability generating functional 

Following Moyal define the probability generating functional (p.g.fl.) with 



X] ~! ( / n^2;^?/'(a;,) I Pn{xi,...,Xn), 

n=l ' \ i=l / 



Gx{i^) =Pa + 2_^—\ I \ \dXi'>p{xi)\pn{xi,...,Xn), (1) 

where Pn(xi, . . . ,a;„) > and 

oo 

PO + ^ I / I Ns^j I P7i{xi, . . .,Xn) = 1 (2) 



£^ ( ^f[dx^]pn{Xl,. 
n=l ' V"^ i=l / 



The p.g.fl. is commonly used in point process theory to characterise point patterns [2]. 
For compactness of notation, we write the p.g.fl. as 



n=0 ' i=l ) 



(3) 



Definition: Gateaux difTerential 

A Gateaux differential |T, p71-74], of a functional ^'('0) with increment f , can be defined 
with 

= + W). (4) 

e-i-O e 

The n*''-order differential (or variation), which is defined by recursively, is denoted 
<5"vi/(v,,a,...,e„). 

Lemma 1 

The p.g.fl. G(K) defines P uniquely, since the Janossy densities can be recovered with 
Gateaux differentials [3], i.e. 

=,5'=G(0; 6, •••,&), (5) 

where ii{x) = ^xAx). 
Proof 

See Moyal [B]. 
Lemma 2 

Let n be the set of all partitions of variables 771 ,...,?]„ , and tt G 11 denote a single 
partition that has constituent blocks a; e tt of size |w| consisting of constituent elements 
Ctt.o;,!, . . . , Ctt.wJwI G {?7i, . . . , ?]„}. The ri*''-order variation of composition fog with 
increments 771 , . . . , r/„ is given by 

-5" ((/ o g) (y) ; Tyi, . . . , 7y„) = ^ (^(y); ^^,1(2/), . . . , [y]) , (6) 

where £,Tr,uj{y) is the variation of order \uj\ with increments C7r,ci;,i, . . . , Ctt,lj,\ui\, i.e. 

■?7r,c^(2/) = '^'"l.g (y; C7r,tu,ij • • • ) ■ (7) 
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Proof 

See Clark [T]. 

Lemma 3 

Leibniz' rule for functional derivatives the differentials of the product f{y)g{y) with 
increments S = {771, . . . , 77„} is 

^" if{y)9{y);m,---,Vn) = '^'*'/(y;'/'i'---''^i*i)'^'"'~'*'5(2/;6,---,'?|H|-i<E.|) , 

(8) 

where the summation is over all subsets $ of the increments, and $ — . . . , 0|$|}, 

s-* = to,---,^iHi-i*|}. 

Proof 

See Mahler p389]. 

2. Deconvolution of point processes 

This section presents the main result. Lemma 4 establishes the higher-order differ- 
ential of the reciprocal of a functional, which is then used to establish a formula for 
the higher-order quotient rule for Gateaux differentials in Theorem 1. It is then shown 
in that the result is directly applicable to point processes. 

Lemma 4 

Consider a functional g{y) in Banach space x- Then the n*^-order Gateaux differen- 
tial [3] of l/g{y) with set of increments 11 = {ryi, . . . ,rin} is found with 

where the summation is over all partitions tt of increment set 11. The product is over 
all cells CO = {Ctt.lj.i, ■ ■ ■ , Ctt.lj.Ilji} of partition tt. 

Proof 

Suppose that f{x) = where x is real- valued and non-zero. Then, following the 
definition of the Gateaux derivative (which in this case reduces to the same rules as 
ordinary derivative), 

<5"(l/x;,7i,...,77„) = i^i^. (10) 
Substituting this into ([6]) gives the required result. 

Theorem 1 

The higher-order quotient rule for Gateaux differentials with increments S = {rji, . . . , ?7„} 
is determined with 
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''{wy"'--")' 

T. (T. ^^S? n (»' ■ ■ .C^-.l-l)) f =1-1"'/ *H-|„|) , 

ncH \7ren wgtt / 

where the summation on the left is over all subsets 11 of S. The inner summation is 
over all partitions tt of subset 11. The term on the right, J'^'g (y; Cir.i^.ii ■ • • ' Cv,u,\u\) is 
the variation of order |E;| — |n| with increments <j)i, . . . , ^|H|-|n| that comprise S — 11. 

Proof 

The result follows from application of Lemmas 3 and 4. 
Corollary 

Consider a point process defined by the superposition of point processes with p.g.fl.s, 
i.e. 

G[h] = Gi[h]G2[h], (12) 

defined with Janossy densities p(xi, . . . , x„), for G[h], g(xi, . . . , x„), for Gi[h], and 
r{xi, . . . ,x„), for G2[h], for n > 0. Let S = . . . ,a;„}. Then we can recover the 
Janossy densities r{xi, . . . , Xn) of G2[h] with 

r(xi, . . . ,x„) = (13) 

Yl i Yl ^ n « K.'^'i' • • • ' Vc.,|c.|) ) P . • . , w'lHi-ini) , (14) 

where the summation on the left is over all subsets IT of S. The inner summation is 
over all partitions tt of subset H. The term on the right p (lui, . . . , u)|H|-|n|) is the 
Janossy density evaluated at . . . , W|s|_|n| which comprise subset S — E, and go is 
the zeroth term in Gi[/i]. When n = 0, ro = Po/<Zo- 

Proof 

Let Gi[ft.] = /(y) and G2\)^ = g(y) in (i). Using Lemma 1, and setting rji = 
Sxi, - ■ ■ ,Vn = ^XnJ and h = yields the result for n > 0. For n = 0, we note that 
G[0] = po, and Gi[0] = go, hence ro = Po/qo- 
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